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We investigate the model dependence of no-helicity flip generalized parton distri- 
bution of the pion upon different approaches for the quark-hadron and quark-photon 
vertexes, in the spacelike region. In order to obtain information on contributions from 
both the valence and non-valence regions, we compare results for spacelike momen- 
tum transfers obtained from i) an analytic covariant model with a bare quark-photon 
vertex, ii) a Light-front approach with a quark-photon vertex dressed through a 
microscopic vector-meson model and iii) a Light-front approach based on the Rela- 
tivistic Hamiltonian Dynamics. Our comparisons lead to infer the same dynamical 
mechanism, the one-gluon-exchange dominance at short distances, as a source of 
both the electromagnetic form factor at large momentum transfer and the parton 
distribution close to the end-points. The expected collinear behavior of the gener- 
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alized parton distributions at high momentum transfer, i.e. a maximum for x ~ 1, 
is also illustrated, independently of the different approaches. Finally a comparison 
with recent Lattice calculations of the gravitational form factors is presented. 

PACS numbers: 12.39.Ki, 14.40.Aq, 13.40.-f, U.lO.St 

I. INTRODUCTION 



In recent years a growing interest in the study of the Deeply Virtual Compton Scattering 
(DVCS) has motivated an impressive amount of work aimed at the extraction of the so- 
called Generalized Parton Distributions (GPD's) from experimental data (see, e.g., Refs. 

for recent reviews). In principle, GPD's allow one to achieve an unprecedented 
level of detail on the knowledge of hadronic states. 

Naturally, the pion GPD should represent a test ground of any approach that addresses 

' hadron structure, and this explains the wealth 
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the issue of obtaining a detailed description o: 

of papers devoted to such a task (see, e.g., llSlp. In what follows, 

we focus on the GPD's that do not depend upon the helicities of the constituents, namely 
we analyze the pion isoscalar and isovector GPD's, as defined, e.g. in 

Aim of our paper is the investigation of the model dependence of those no-helicity flip 
(chiral-even) GPD of the pion upon different relativistic approaches in the spacelike region, 
i.e. for negative values oft= {p' — p)^, where p and p' are the initial and final four-momenta 
of the pion, respectively. In particular, the study of the GPD's in the valence and non- 
valence regions (see the following Section) is emphasized by the choice of three different 
models that explore different kinematical regions: i) a covariant analytic constituent quark 
(CQ) model, that covers the whole kinematical domain and allows us to interpolate between 
the other two models; ii) two phenomenological models, elaborated within a Light-front (LF) 
framework (see e.g., [la, [l7|, [iSj for a review), which have a smaller kinematical range of 
applicability, namely one addresses the non-valence region and the other the valence one. 

The first model is analytic and covariant, and depends upon the mass of the constituents 
and a parameter, fixed by the decay constant of the pion. The main ingredients of such 
an approach are: i) the Bethe-Salpeter amplitude (BSA) of the pion, modeled through 
an analytic Ansatz in the Minkowski space, ii) the Mandelstam formula 191] (or Impulse 
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Approximation formula) for the matrix elements of the current operator, and iii) a bare 
quark-photon vertex. A peculiar feature of our Ansatz for the pion BSA is given by the 
symmetry under the exchange of the constituent momenta. A first version of such a model 
was adopted in Ref. 20|] to investigate the frame dependence of the description of the 
electromagnetic (em) pion form factor, putting in evidence the possibility to study the non- 
valence content of the pion by using a suitable reference frame. In the present work, we 
consider a natural extension of the model, that features a better end-point behavior of the 
BSA, as well. 

A second model, developed within the LF Dynamics and already applied to the em pion 
form factor in both the space- and timelike regions 2lj], is still based on the Mandelstam 
formula. However, this model retains only the analytic structure given by the poles of the 
Dirac propagators in the analytic integration over k~ = — k^, i.e. the minus component of 
the constituent four-momentum appearing in the loop formula. An important consequence 
of the ^"-integration can be reached in a frame where the plus component of the virtual- 
photon four-momentum is different from zero, i.e. A+ = A'^ + A'^ 7^ 0. Indeed, in this 
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lowing 



23| for 



frame the contributions in the valence and non-valence regions can 
an investigation of the Fock components of the hadronic state (see 
an overview of the Fock expansion of a hadron state, within the LF framework). Another 
relevant feature of this model, that has a fundamental impact in the timelike region, is the 
quark-photon vertex dressed by a microscopic version of the vector meson model (VMD) 
2l| . Finally, as explained in detail in 21[ , the model lives in the non- valence region, in the 
limit of vasnishing pion. 

A third model is constructed within the LF Relativistic Hamiltonian Dynamics (LFHD), 
where the Poincare covariance is fully satisfied (see, e.g. for a detailed review). In par- 
ticular the rotational covariance is fulfilled through the introduction of the Melosh rotations 
and the proper definition of the total intrinsic angular momentum. At the present stage, 
the model explores only the valence region. 

The paper is organized as follows: in Sec. [Ill a brief resume of the general properties 
of the pion isospin-dependent GPD's is presented; in Sec. Illlt the Fock decomposition of 
the GPD's is discussed, in view of a frame-dependent analysis; in Sec. IIVI a covariant CQ 
model, that allows an analytic evaluation of the pion GPD's is described; in Sec. IV A[ a 
first CQ Light-front model, with a quark-photon vertex dressed by a microscopic version of 
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the vector meson dominance model, is presented; in Sec. IV Bl the LFHD model, where the 
full Poincare covariance is implemented is described. Finally in Sec. IVII and Sec. IVIII the 
results are discussed and the conclusions drawn. 

II. PION GPD'S: KINEMATICS AND GENERAL FORMALISM 

In the spacelike region, let us first illustrate the kinematics of the DVCS process with 



the symmetric momenta convention shown in Fig. [T] (see 12| for the reduction of the DVCS 
diagram to the one presented in Fig. [H and the pioneering paper 2J| for the DIS regime), 
i^'or on- mass-shell pions, i.e. p'^ = = m^, and adopting standard notations (see, e.g. 

13) 

e = -^^ = -^ = ^^^ m<i) 

^ 2 P-n 2P+ p++p'+' vi',1 - 

x = -^ = -^, (l>x>-l) , (1) 

where n is a light-like 4- vector, such that = n ■ v = v*^ + (the scalar product is defined 
as a ■ 6 = {a^h~ + a~6+) /2 — • b^), P = \{p' +p), and k is the average momentum of the 
active quark, i.e. the one that interacts with the photon (see Fig. [T]). Notice that p"*" and 
p'~^ are necessarily positive, while A"'' > is taken by choice. From Eq. ([T]) one can trivially 
obtain the following useful relations 

As it is well known, the variable x allows one to single out i) the valence region (where one 
has only contributions diagonal in the Fock space , cf the following Sec. IIIII) given by the 
union of two intervals: x G [—1, — 1^|] (corresponding to an active antiquark) and x G [|.^|, 1] 
(corresponding to an active quark), and ii) the non-valence region, x G [— |^|]. In Fig. 
[2]^a), it is shown a representative of the contribution with an active quark in the kinematical 
region x G (all the constituents have a plus-component of their-own momentum 

bounded from above by the corresponding quantity of the parent pion). In Fig. EJ^b), it is 
shown a contribution from a pair-production process, non diagonal in the Fock space. In 
Appendix |A] a more detailed kinematical discussion is given. Finally, as a short detour, let us 
remind that the pion BSA, integrated over the minus component of the quark momentum. 
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yields the two-body Fock contribution to the pion state, notably non vanishing only in the 



valence sector (see 



18|) 




p' = P + ^ 



k-P 



p = P 



A 
2 



FIG. 1: Diagrammatic representation of the pion GPD, with four-momenta definitions. 



p' = F + f 





FIG. 2: LF time-ordered analysis of the pion GPD. Diagram (a): a contribution in the valence 
region, 1 > x > |^|, (see text). Diagram (b): a contribution in the non-valence region, |^| > j; > 
— (see, text). Vertical dashed line indicate a given value of the LF time, in order to single out 
the number of constituents in flight. 



ence region is called DGLAP [25|] region, 
I27I region. 



Within the QCD-evolution framework, the va 
while the non-valence one is called the ERBL [26 

In the interval [|^|, 1], the relation between the LF momentum fraction, Xq, of the active 
constituent in the initial pion (with the support [0, 1]) and the variable x defined in Eq. ([T]), 
is given by 

k+-A+/2 A;+-A+/2 x + ^ x - \^\ 



(3) 



The isospin-dependent GPD's (see, e.g. [4, [tI, [ill, [l3|]) are the matrix elements of light 



cone bilocal operators separated by a light-like distance, z"^ 



, evaluated 
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between pion states with different initial and final momenta. In the light-cone gauge, where 
^giuon -^1 = and the gauge link becomes unity, one can introduce isoscalar and isovector 
combinations for the off- forward (t 7^ 0), non-diagonal 7^ 0) GPD's, as follows 



dz 
Air 



JxP^ 



5=0 



(4) 



and 



nP+z 



z=0 



where z = {z~^, z±}, while i'qiz) and r^ipq^z) are the following doublets of quark field 



(5) 



respectively. In Eqs. 



J4]L and 





(6) 

following [3], instead of the Cartesian components, 
7r°, TT^, TT^, (adopted in [4, [l3|), the charged pions have been introduced, viz 



k ) = k ) • 



(7) 



The functions H'^ix, ^, t) and H'^{x, ^, t) are u and c? GPD's, respectively, and contain quark 
and antiquark contributions (cf. the parton interpretation of if, e.g., in l|, |J] and Fig. 2). 
It is worth noting that H'^{x,^,t) has the support x G [—1, 1]. Finally, due to the isospin 
symmetry one has 



TJU Tjd 



and combining charge and isospin symmetry (G-parity) one gets 

H:^ix,^,t) = -H:.i-x,^,t) = -H^^{-x,^,t) . 



(8) 



(9) 



In what follows, we deal with a charged pion and the subscript vr"*" in the quark GPD's is 
dropped out whenever no ambiguity is present. 

For vanishing 1^ and t, one has the following partonic decomposition (cf. [l. 



H'^ix, 0, 0) = 9{x) u{x) - 9{-x) u{-x) 
H\x, 0, 0) = 9{x) d{x) - e{-x) d{-x) 



(10) 
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Equations ([8]) and (Q together with the partonic interpretation lead to the well known 
relations between the standard parton distribution functions (let us remind that the relations 
pertain to active quarks), viz 

Ut,+ {x) = d^-{x), (x) = (J^+ (x) . (11) 



The symmetry property of H^^^'^{x,^,t) (see, e.g. 0, Q]) under the transformation 
X —X, that just reflects i) the charge-conjugation {p — > —p and p' — > —p') and ii) the 
isospin invariance, reads (reminding Eqs. ([8]) and IQ) 

H'=\x, i,t) = ]^ [H\x, e, t) - H^i-x, t)] = -H'=\-x, e, t), (12) 
H'=\x, t) = ^ [i^"(x, e, t) + H\-x, t)] = H'=\-x, e, t)) . (13) 

Therefore the two GPD's are odd or even in x depending upon the isospin combination. 
In addition, under the transformation ^ —* — ^, that amounts to apply a time-reversal 
transformation (since we have to exchange the initial and final pion momenta) and to exploit 

riQ 

Hermiticity, one has (see, e.g. [4,171]) 

H\x,^,t) = H'{x,-U), (14) 

namely H^{x,C,,t) must be even in ^. 
From Eqs. (jl]) and ([5]) one has 

H^x, e, t) = H'='{x, e, t) + H'=\x, e, t), 

H\x,i,t) = H'=\x,i,t)-H'=\x,i,t) . (15) 

As well known, the following sum rules hold (note a different overall factor with respect 
to [nl, Q due to our choice of dealing with a charged pion, cf Eq. ^^^) 

J dxH^=\x,^,t) = j dxH''{x,^,t) = F^{t), (16) 

j dx X H'=\x,^,t) = j dx X H'^ix.i.t) = ^ [^2(t) -e'^i(t)] . (17) 

n Eq. ( [T6l) . F.„{t) is the pion em form factor (see Appendix [C]) , while, according to the Ref. 
28], in Eq. (11 7p Oi{t) and 62{t) are the gravitational form factors (see also, e.g.. 



that enter in the parametrization of the matrix elements of the quark part of the energy- 
momentum tensor (notice that in the chiral limit one has 6'i(0) —^2(0) = (9(m^)). It should 
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be pointed out that the sum rule (fTTIl for t = and ^ = yields the longitudinal-momentum 
sum rule for the pion, i.e. < Xg >, as numerically illustrated in Sec. IVll 

For vanishing ^ and t, one can exploit i) Eqs. (fT2i) and (|T3l) and ii) the partonic decom- 
position (cf. Eq. ( ITTl) ) obtaining 

H'=%x, 0, 0) = i [H^{x, 0, 0) + H'^ix, 0, 0)] = ^ [H^{x, 0, 0) - H^i-x, 0, 0)] 

= 9{x) ^ [u{x) + u{x)] - 9{-x) ^ [u{-x) + u(-a;)] , (18) 



and 



H'=\x, 0, 0) = ^ [H'^ix, 0, 0) - 0, 0)] = i 0, 0) + 0, 0)] 

= 9{x) - [u{x) -u{x)]-0{-x) - [u{-x) -u{-x)] . (19) 
2 2 



Analogous relations, with singlet, q{x) + q{x), and valence, q{x) —^{x), combinations, for 
the (i-quark can be easily obtained, by using Eq. ffTTl) (see also [7|). It is worth noting 
that for ^ = A+ = the ERBL region shrinks to zero and the variable x reduces to Xq 
(Eq. (j3])). Finally, from Eq. (fT6|) one has a normalization for the valence combination 
Uv{x) = u{x) — u{x) given by dx u^{x) = 1. 

It should be pointed out that the parton distributions represent a bridge toward the chiraL 



even transverse- momentum dependent (TMD) distribution, /i(x, |kxp) (see, e.g. |29l. I30l.l31| 
for the nucleon case), as shown by the following relation 

q{x) = J dk^ /f(x, |k^n, (x > 0) . (20) 

Furthermore, it is worth noting that an experimental access to fi{x, |k^p) and to other 
TMD's is a fundamental step in order to understand the correlations between constituents 
inside the pion, and eventually the dynamics. 

To complete this brief resume of the general formalism, we have to mention that the 
sum rules in Eqs. (1161) and ( |T7I) . are the lowest order of the moments of the isovector and 
isoscalar GPD's. In particular, H^^^{x,C,,t) (see Eq. fll3p ) has only even moments, while 
H^^^{x, ^, t) (see Eq. ( IT2!) ) has only odd moments. Moreover, it turns out (see, e.g., [4|) that 
the n-th Mellin moments of the GPD's are polynomials of C, with highest power n for even 
moments and n + 1 for odd moments, i.e. only even powers of ^ appear, as expected from Eq. 
( IT4l) . It is worth noting that the so-called polynomiality follows from general properties, like 
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Hermiticity, covariance, parity and time-reversal invariance The isospin-dependent 

moments are given by (j > 0) 



r dx x'^ H^=\x,^,t) = i2^wM20'\ 
/ dx x^^^' H'=\x, ^,t) = J2 Ai=^2,2^{tmf ■ 

J -I ,-n 



(21) 
(22) 



In particular, numerical calculations of i) F.,^{t) = A{Q^{t) and ii) A2q^ = 92(t)/2 and 



i/=o 



A^=o = _^^(t)/8, will be presented in Sec. |VT1 

In conclusion, approaches that satisfy the basic field-theoretic assumptions underlying 
polynomiality, i.e. extended Poincare covariance, automatically fulfill the conditions (12T|) 
and (l22l) . In general, such a property is an important test of consistency of the model. 

III. FOCK DECOMPOSITION 



Let us introduce the Fock expansion of the pion state, taking care of the colorless feature 



of each components and inc 



in this Section, (see, e.g. jl7,Q), viz 



uding the amplitudes inside the kets to simplify the notations 



k) = \qq) + \qq;g) + \qq; qq) + ... 



(23) 



Then one can decompose the GPD's in terms of their Fock contents (see also [4, 122]), i.e. 
one can rewrite Eqs. (jl]) and ([5]) by using e.g. 

Fock 



i7^(a:,e,t) = Y,{n;n\V%\'K-n)e{\x\ - \i\)e{l - \x\) + 



+om-\A) 



n 

'Fock 



Fock 



^(vr; nlr^^lvr; n + 2)^(0 + J](7r; n + 2|r^^|7r; n)e{-0 



+ ... ,(24) 



where n indicates the number of quarks and antiquarks, Fd and Fnd are the diagonal and 
non diagonal, in the Fock space, terms of the current operator, and dots represent all the 
other transition matrix elements, possibly containing states with gluons. The diagonal terms 
yield contributions to the valence region (DGLAP region), while the non diagonal ones have 
to be considered in the non- valence region (ERBL region). In Eq. (l24ll . we have shown only 
transitions involving fermionic fields, and this explains the selection rule An = 0, 2. 
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In a simple picture of a hadron, the valence state has a dominant role at the hadron scale, 
and this leads naturally to associate the DGLAP region with this Fock component. 

The same decomposition can be applied to the em and gravitational form factors, and to 
all the t-dependent "generalized" form factors appearing in Eqs. ( 1211) and ( 122|) . Clearly, this 
kind of decomposition could allow a deeper understanding of the dynamics related to the 
components beyond the valence one. As a simple application, let us consider the em form 
factor. From Eqs. f|T6|) and fl2^ . retaining only the fermionic transitions, one has 

pl Fock pi 

Fi^\^,t) = 2 dxH'=\x,i,t)=2y^ dx (tt; nir^rV; n), (25) 
r\i\ 

F(™)(e,t) = 2/ dxH'=\x,i,t) 
Jo 



Fock 



= 2V / dx [6(0 (vr; n|r^=^V; n + 2)+e{-0 (vr; n + 2|r^7,V; ^)] • (26) 

The valence term, Fn'"\^,t), receives the largest contribution from the valence component 
of the pion state, but it does not give the full result in the whole kinematical range, as 
indicated by the residual dependence upon ^. The non- valence term, Fi-"^^(^,t), is due to 
contributions like the pair-production mechanism, see Fig. Mjo). The sum of Eqs. (|25ll and 
(!26|) leads to the full result, viz 

F^t) = Fi^\U) + Fi^^\U) (27) 

and it is independent of ^ and function of t only. One can also express the invariance of the 
sum under changes of ^ as: 

am am 

;Fi''\U) = -^Fi-^\U) , (28) 



Q^m TT v^w 

with m > 1. It is worth noting that all the derivatives of F.„^{t) are independent upon ^, and 
therefore relations like the one in Eq. fl28|) can be generalized, i.e. 



am a£ am a£ 

with m > 1 and £ > 0. As a consequence, with the help of Eq. ( |T6l) . one can deduce 
interesting sum rules for the partial derivatives of H^^^{x,C,,t). 

Let us remind that calculations of the elastic form factors have been performed in different 
frames. In particular, it has been chosen i) the Drell-Yan frame, where A"*" = and therefore 
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^ = (see Ref. [17|] for generalities on the Drell-Yan frame), or ii) a Breit frame (i.e. 



= —A ) where A i =0 (see 



32] for an extended discussion of the motivations for 



adopting such a frame), and then follows a kinematical trajectory in the {C,, t)-plane given 



by 1^1 = l/yl — 4:m'^/t (see below Eq. f l55]) ). In the first case, the em form factor is 
saturated by the valence contribution, because of ^ = (cf Eqs. (|27j) and (1251) ). while in the 
second frame both valence and non-valence terms contribute, since ^ does not vanish, but 
changes with t. For — >> the value of C, approaches 1, and therefore the non- valence 
term saturates the em form factor (cf Eqs. (1271) and fl26|) ). In model calculations this general 



behavior was indeed observed 



20l ]. It is understood that for an experimental investigation of 



the whole (^, t)-plane, different kinematical conditions are needed, also exploiting the helpful 



properties of the LP boosts (see, e.g. jig). 

Following the same spirit, one could extend this analysis to the other form factors that 
appear in Eqs. (12T!) and fl22l) . i.e. one can consider the partial derivatives of the valence and 
non- valence contribution to the generalized form factors ^2j+2,2i(^) ^2j+i,2i(^)) obtaining 
final relations that have the same structure as the ones in Eqs. (1281) and (1291) . 



IV. COVARIANT MODEL OF THE PION WITH PAULI-VILLARS 

REGULATORS 



In Ref. [20|, an analytic covariant model, symmetric in the exchange of the constituent 



four-momenta (see Refs. 



33 



3J] for previous non symmetric versions) was adopted for 



evaluating the em form factor. In this work, a direct extension of the symmetric covariant 
model to DVCS is exploited for calculating the no-helicity flip GPD, in the spacelike interval 
> t > -10 (GeV/c)2. 

In a Breit frame, one has A° = 0, i.e. A+ = — A^, and p' = — p = A/2. By choosing 
A"*" > 0, and reminding that 

,„ m2 + |Ax|V4 _ m2 + |A^|V4 
P = J- , P = , 

A- = p'~ -p- = -A+ " ' ^' ^ = A- ,\ \t\ L ^ 30 

p'+p+ (p+ + A+)p+ ^ ^ 



one gets 



-A+ + J-A^ + Ami l^ 

^ = -2 -= - -Y^'^-C 
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,^ A+ + y/- + Ami A+ 1 



Then the following relation holds (notice that 2 P+ = -y/^AM-^mf ) 

A2 = _A+2 _ I A^|2 = -4e'p+2 _ I A^|2 = (-A2 + Ami) - (32) 

which leads to a constraint on the maximal value for the variable ^. As a matter of fact, in 
the spacelike region — A^ + Am^ ^ and one has 

_ -A^-|A,P 

Then, the maximum value of is found for = 0, viz 

-A^ 

< — = — < 1 . (34) 

^ - -A2 + Ami - ^ ' 



For m,r = (and A 7^ 0), one has 



|A 



|2 



e^ = l + ^- (35) 



If one additionally chooses a frame where = (i.e. only A^ 7^ 0), then ^ = — 1 and 
therefore, in this extreme case, only the non-valence region contributes. 

A basic ingredient in the analytic covariant model of Ref. j20j is the pion BS amplitude, 
that can be quite well approximated by retaining only the pseudo-scalar Dirac structure 



[see, e.g.. 



35|). namely 



m{k -P,p) = -— Sik- A/2) 7^ A{k - P,p) Sik- P), (36) 

J-rr 

wher e ^ is the quark-pion coupling, as suggested by a simple effective Lagrangian (see, 
e.g. (2J]), /tt = 92.4 MeV the pion decay constant, m and S{k) are the mass and the 
Dirac propagator of the constituent quark (CQ), respectively. In Eq. (|36ll . A{k — P,p) is a 
scalar function that describes the momentum-dependent part of the coupling between the 
constituents and the spin-0 system and plays the role of the Pauli-Villars regulator of the 
otherwise divergent integrals that yield GPD's or the em form factor. In particular in this 
work we adopt two symmetric (in the exchange of the CQ four-momenta) covariant forms: 



i) the one considered in Ref. 201], and based on the following sum 



Mk-P,p)= C, {t7—-—^ ^+ r,p .J 2 ^ 1 ^ (37) 

[(A; — A/2) — m|j + «ej \[P — k) —mj^ + ie\ 
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and ii) a natural extension based on a product, viz 

1 1 

AafA; - P,p) = Cs 7 5 t t 9 t ■ (38) 

This product-form provides a more realistic transverse-momentum fall-off, as seen from the 
expected behavior of the BS amplitude obtained by using a simple (one-boson-exchange) 



kernel (see, e.g., [36||)), and this has a sizable impact on both the high- momentum tail of 
the em form factor and the end-point behavior of the parton distribution, as shown in the 
results presented in Sec. |Vll We can anticipate that the most favorable comparison with 
the experimental data of the em form factor is obtained by using the product-form, as also 
expected if one follows a pQCD analysis, where a one-gluon exchange represents the leading 



contribution to the kernel 



'27 



37|. 



In both expressions, once the constituent mass m is chosen, m^i is determined by fitting 



the experimental value for f^^ (cf j20|), while the constants Ci and C2 are fixed by exploiting 
the charge normalization, as discussed below. 

As a final comment on the Dirac structure that appears in Eq. fl36l) . we remind that it 
leads to the standard Melosh rotation for a pair of fermions coupled to a total spin = 



(see |38|), once we consider the valence wave function, defined as follows (see, e.g. [l8|) 

^val{l^'^ ,l^±,P) = [ Son{fi-p) 7^ A(k,p) Soni-H) (39) 

In J ^TT 

where S'on(/t) = (^on + rn)/ {k^ - m^ + it) with = = (m^ + \k±\^)/k^, k+, 

The no-helicity flip GPD for the pion is calculated in one-loop approximation (triangle 
diagram cf. Fig. [1]) with the BS amplitude of Eq. (136|) and the symmetrical forms shown in 
Eqs. ( 1371) and ( l38i) . In particular, the w-quark GPD is given in Impulse Approximation by 

X J ^^AP'''^ - V^{k,p,p')A{k - P,p') A{k -P,p), (40) 

where A^^^ = 3 is the number of colors, TZ = 2m^//^ and 

V^{k,p,p')=Tr!^S{k-P)^' s(^k + ^^^+ s(^k-^^^'j . (41) 

The presence of the delta- function in Eq. fl4Up . given the kinematical relations in Eq. ([T]), 
imposes the correct support [— 1] for the variable x as discussed in details in Appendix IB] 
(note that H'^ has the support [—1, |,^|] for the variable x) . A relevant feature in the analysis 
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of the GPD, as well as in the case of the em form factor, is given by the instantaneous term 
present in S{k). As a matter of fact, the Dirac propagator can be decomposed using the LF 
kinematics as follows 

'^^^^ - p _ ^2 + - ^on{k) + 2k+ - k+ik- - fc- + ^) + 2A;+ ' ^^^^ 

where the second term, proportional to 7^, is an instantaneous one in the LF time. It 
should be pointed out that the instantaneous contribution to the GPD is produced only by 
the spectator fermion (in the present example an antifermion), i.e. by S{k — P). Indeed, the 
instantaneous terms pertaining to the other propagators do not contribute, because of the 
property (7"*")^ = 0. In our symmetric model, the instantaneous term of Eq. (142!) contributes 
to H^{x,^,t) both in the valence and in the non-valence region (see Eqs. ( 1B10I) -( |B15[1 ). 
since we take fully into account the analytic structure of the symmetric vertex function (for 
a different approach, where such an analytic structure is disregarded see 
The pion em form factor is obtained by using the sum rule f|T6|) : 

= -^^c ^p>-Jlp+^ I V^{k,p,p')Aik-P,p') A{k-P,p) . (43) 

The last expression for -F^(t) can be extracted directly from the Mandelstam formula for 
the matrix elements of the em current [3] (see, e.g., 3^, 0]); well. Notice that the 
model preserves current conservation, as discussed in 201 ]. 

The normalization of the form factor, Eq. (1431) . allows us to determine Ci and C2 in Eqs. 
( !37|l and (138|) . Such a charge normalization represents the irapulse approximation of the 
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4l|. 



normalization condition in the fully interacting BS theory 

A standard analytic integration on k~ (see Appendix[B]for details) leads to the following 
decomposition of H^[x,C,,t) in valence and non-valence contributions 

H-{x,C,t) = H}^,){x,C,t)9{x- 1^1)9(1 -x) + H};^,^{x,^,t)9{\^\-x)em+x) . (44) 

Notice that if^" ^ and H^^^-^ are given in Appendix [B] for the two momentum dependences 
shown in Eqs. f l371) and fl38l) . 

The (i-quark GPD can be obtained reminding Eq. Q. 

Within our covariant model the valence component if^^-j in Eq. fl44p is an approximation 
to the diagonal terms in Eq. (1241) . while the component H^^^-^ contains the contribution 
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of the pair-production mechanism from an incoming virtual photon with A"*" > and 
approximates the non diagonal terms. 

An interesting approximation of the contribution to GPD in the valence region can be 
obtained once the analytic structure of the BS amplitude is disregarded and only the poles 
of the propagators are retained in the integration over k~ (see Appendix [B]) . As a matter 
of fact, see Eq. (]B7p . within the mentioned approximation 



K+{p+ - K+){p'+ - K+) 



where k = P — k, = {w? + and 

Tr[0^{^^-n)] = + m) [(/ - + m] 7+ [(^ - + m]} . (46) 

Moreover, if in Eq. fHSl) we identify the following ratio 



with a model LF wave function, then the final expression coincides with the result obtained 

owing Sec IV Bp . since the trace Tr[0^{K~^)] generates 



within a LFHD approach (see the 
the correct Melosh-rotation factor 



bi: 



. We would stress that the identification is meaningful 



once the analytic structure of the BS amplitude is disregarded. 

V. LIGHT-FRONT MODELS OF THE PION 

In this Section we present models that at different extent i) fulfill the Poincare covariance 
and ii) take into account the Fock components of the pion state beyond the valence contri- 
bution. A first important difference between the models is given by the frame we choose. In 
the approach we call Mandelstam-inspired LF model, a Breit frame, where Aj^ = 0, is con- 



sidered. This choice was followed in Ref. 2l|] in order to perform a microscopical calculation 
of the em pion form factor in both the space- and timelike regions. It should be pointed 
out that such a frame leads to consider contributions from a pair-production mechanism, 
differently from what happens in a Drell-Yan frame, where = 0. This second frame is the 
one adopted in the second approach illustrated in this Section, based on a LF Hamiltonian 

n 

Dynamics description of the pion state (see, e.g. p^] for a general review of LFHD). 
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A. Mandelstam-inspired LF Model 



In Ref. 2l|] an approach was elaborated to calculate the em form factor of the pion 



starting from a covariant expression of the matrix elements of the current given by the 



Mandelstam formula 191] (cf also Eq. (HOl) ). Moreover, a microscopic VMD was used for 
dressing the quark-photon vertex. The dynamical inputs of such an approach were the wave 
functions of both the pion and vector mesons, taken as eigenstates of the relativistic CQ 
square mass operator of Ref. which includes both confinement, through a harmonic 

oscillator potential, and vr — p splitting through a Dirac-delta interaction in the pseudoscalar 
channel. In what follows we apply the same approach for evaluating the no-helicity flip 
GPD's. 

Let us first illustrate the kinematics in the adopted frame, where A_|_ = (i.e. = 
A^/A+) and = p'j_ = 0. Then in the spacelike region, for A'^ > 0, one has for and 




smce 



2 2 

mi mt 



p'+ p-i 



2 

mz 



A- = p'^ -p- = -A+ — . (48) 

{p+ + A+) p+ ^ ' 



The following simple relation between ^ and A^ holds 



^— (49) 



It is easily seen that if vHt^ = one has ^ = — 1 for any A 



2 



Extending the approach of Ref. [2l|, one can find for the quark GPD the same formal 
expression of Eq. fHU]) . but i) a microscopic VMD dressing, V^lk, A), is considered instead of 
the bare quark-photon vertex, 7^^, and ii) phenomenological Ansatzes for the BS amplitudes 
in the valence and non-valence regions are adopted. Another basic difference with respect 
to the analytic model presented in the previous Section, is that only the simple analytic 
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structure of the Dirac propagators is retained, i.e. the analytic structure is disregarded in 
the BS amphtudes of both i) the initial and final pion and ii) the VM dressing of the quark- 
photon vertex. This approximation turns out to be a very effective one in the calculation of 



the em form factor just in the = frame 

In Ref. 2l|], a further simplification in the calculation was achieved by a quite natural 
assumption, namely a vanishing pion mass. Within such an approximation only diagrams 
with a qq production contribute (cf Fig. [2]^b)), and this implies the necessity to introduce 
the VMD dressing. We have to stress that a bare term is missing, due to the vanishing pion 

n 

mass (cf the discussion in [2 11]). Therefore, in the quark-photon vertex for the covariant 
model, Eq. (14T1) . the Dirac matrix 7+ is replaced by the plus component of the following 
four-vector, that microscopically describes a VM dressing. For t < one has 



Vn 



(t - Ml) 



(50) 



where fvn is the decay constant of the n-th VM into a virtual photon {calculated in the 
model), P^ = {Af^/A+, A+, 0_l} the four-momentum of an on-mass-shell VM with a square 
mass given by P^ = and e\{Pn) its polarization. Moreover, the VM BS amplitude is 
approximated as follows 



'i/nx{k,Pn) 



^ + m 



ex{Pn)-Vn{k,Pn) Kn{K Pn) 



Pn + m 



(51) 



{k — Pn)^ — + te 

where Vn{k, Pn) is the proper Dirac structure, and A„(fc, Pn) the momentum-dependent part, 
approximated on the LF hyperplane, as discussed below. 

In the valence sector, after performing the k~ integration, both pion and VM's BS ampli- 
tudes reduce to 3D amplitudes with one constituent on its mass shell. In 2l|, the momentum- 
dependent part of the on-shell VM BS amplitude (that contains on both sides proper Dirac 
projectors) is described through a LF VM wave function, i.e. 

A.(fc,P„) 



[Ml-M^{k+,kr,P;t)] 



V^„(fc+,k^;P+), 



(52) 



and 



Mo2(/c+, k^; P+)] = P+ [k;^ + (P„ - k)-^] . 

In Eq. (1521) ■ ipnik~^, kj,; P+) is an eigenfunction of the relativistic CQ square mass operator 
of Ref. 42|, as mentioned at the beginning of this Section. Moreover, it is normalized to 
the probability of the valence Fock state, according to the model elaborated in [21|. 
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The valence component of the pion was modeled adopting an analogous Ansatz. More- 



over, m 



2l| two different calculations were generated by using i) the pion eigenstate of the 



model in Ref. 421] and ii) the pQCD asymptotic wave function (see, e.g. 27l|). 

In the non- valence region, namely the only region contributing to the GPD's for m,r = 
(see Eq. (H9l) ). besides the pion valence component in the initial state one has to deal with 
a non-valence component of the pion state, since the process depicted in Fig. Mjo) can be 
interpreted as a transition from a state composed by the valence component of the initial 
pion and the virtual photon, \qq, 7*), to a higher Fock component, \qq, qq), pertaining to the 
final pion. At level of the pion BS amplitudes, one has to model an off-shell BS amplitude, 
that takes into account the absorption of the initial pion by an antiquark (according to the 
case illustrated in Fig. [2](b)). In 2l|] a simple Ansatz, namely a constant vertex was assumed, 
like in Ref. Notice that such a coupling constant is fixed by the normalization of the 

pion form factor, since the diagram shown in Fig. [2]^a) does not contribute, as a consequence 
of the simplification = 0. 

Within the approach presented in this subsection, since = 1, (given the vanishing m^) 
the quark GPD has only contribution from H^^^y i.e. 



if"(a;, \i\ = 1, t) = i/("„,)(x, \i\ = 1, t) 9{l - x) d{l + x). 



(53) 



where, introducing k = P — k, 



H^nv){x, 1, t) ^ ^ _ ^2 (2vr)3 72 ip+ {p'+ - K+) ip+ - K+) 



where is the constant describing the off-shell quark-pion vertex, while Xi and I2 are given 

by 



1 m 

2 U 



A((p'-k),p')L-=p'— 



X Tr <^ 7+ [(/ - fi)on + m] Vnz{p' - P„) [{^ - ^)o„ + m 



X2 = -Tr I (^o„ + m) [(/ - ^)on + m] Vnz{p' - Pn 



7 



(55) 
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The Dirac structure, 



where all the constituents are on their own mass- 

Fur- 



shell, is chosen in order to generate the proper Melosh rotations for ^Si states 
thermore, the traces previously shown contain the instantaneous terms (see Eq. (142|) ) that 
survive after assuming m,r = 0. In order to model the instantaneous part of the vertex 
functions directly attached to 7"*", we performed the following replacements 



m 



-A((p'-«:),p')L- 

Jit 



-{p'-K,)on 



(56) 



p 



for the pion, and 



A„((p'-/s:),P„)|,_ 



=p' -{p'-K)on 



CvM l/^niip' - I^V, -l^±; Pn 



(57) 



for the VM's, as in [2l|. In Eqs. (!56|) and (l57j) . the constants and Cvm roughly describes 
the effects of the short-range interaction. Indeed, a relative weight, wvm = Cvm/C-k, can be 
used as a free parameter. Let us remind that the on-shell part of the BS amplitudes have 
on the left and right sides the proper Dirac projectors. 

Finally, it is worth noting that the results presented in the following Section IVII have 
been osculated by u^ing alUhe parameters adopted in fl, bnt wtth a CQ ma. » = 200 
MeV and wvm = — 1 (see [211] for m = 265 MeV and different values for wvm)- It should be 
pointed out that only one adjusted parameter is necessary for describing the em form factor 
in the spacelike region. 

The model remains invariant for kinematical transformation, after the approximation we 
have apphed. 



B. Light-front Hamiltonian Dynamics model 



introduced by Dirac in 



Within a LFHD approach (see [16|] for a review of the three forms of the relativistic HD 



45] ) the Poincare covariance of the description of the pion can be fully 



implemented, once the current operator is chosen in order to fulfill the proper commutation 
rules with respect to all the generators (i.e. both the kinematical and the dynamical ones). 
A widely adopted strategy, within the LFHD approach, is to model the em current by using 
a one-body operator, but in the Drell-Yan frame, i.e. where = 0. For instance, in 
this frame the em form factor can be obtained by using only the matrix elements of the 
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plus component of the current operator, and this allows to overcome some difficulties that 



and 



32| for a general 



manifestly appear for hadrons with angular momentum > 1 (see 
discussion). 

In the Drell-Yan frame, A_l 7^ and one can choose p'j_ = —p± = A^/2. It is worth 
noting that only the spacelike region can be addressed, since = — |A_Lp. Moreover, one 
has = p'^ and therefore ^ = for any A^. 

In this section, the LFHD model with CQ's, alread y su ccessfully applied for describing 
the charge form factor and decay constant of the pion 



46 



47j, is adopted for investigating 



the DGLAP contribution to the no-helicity flip GPD. This corresponds to consider in the 
Fock-space expansion of Eq. (1241) the diagonal contribution with n = 2 constituents (i.e. 
the valence component, cf also Eq. ( 139|) . introducing the explicit representation in terms 
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22l |. The quark contribution to the 



of overlap of light-cone wave functions (LCWFs) 
GPD in the region < x < 1 can be written in terms of the LCWF \l/7r(x, k±; Ag, A,) for the 
quark-antiquark system as 

dk I 



X 



u{x, k_L + A' ) u{x, kj 



k+ + 



A+ 



7 



k+ - 



A+ 



^^(x, k;_l; Ag, \i 



E 

Ag , Ag 



dK,± 
2(2^ 



%{x, Ag, Ag) ^'^(x, Ag, A 



9/' 



(58) 



where u(x, k;j_. A) is a LF Dirac spinor (see, e.g. |38|), and Aj are the spin projections. The 
perpendicular component of the active quark momenta, ± Ax/2, become in the intrinsic 
frame 



'1 



X 



kx + (l 



X 



+ (1 - x) A 



(59) 



2 ' ^ ^ V . 2 

with X given by Eq. ([1]). Notice that in the Drell-Yan frame Xg = x (cf Eq. ([3])), since ^ = 0. 

For the model calculation, we use a phenomenological LCWF which satisfies Poincare 
covariance and is eigenstate of the total angular momentum operator in the Light-front 



dynamics. As outlined in Ref. 



461, these properties can be fulfilled by constructing the 



wave function as the product of a momentum wave function '?/'(x, i^±), which is spherically 
symmetric and invariant under permutations, and a spin wave function, which is uniquely 
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determined by symmetry requirements. Therefore, within LFHD one has 
where k = Hz} with 



and the free mass defined by 



= Mo(x, (x - (61) 



The spin-dependent part contains the Melosh rotations Rm{i^) which convert the instant- 
form spins of both quark and antiquark into LF spins and ensure the rotational invariance 
of the pion wave function. The representation of the Melosh rotation is explicitly given by 



For the momentum- dependent part of the pion wave function we adopt the following 



exponential form used in Refs. 
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47| 

1/2 



^v(..«x) = |2(2.)T'^(^^|f^) ;pi7^-p(-.V(2/3^)). (64) 



The wave function in Eq. flMl) is normalized as 

(reminding that dnz = dx Mo{x, «;_l)/[4x(1 —x)]), and depends on the free parameter (3 and 

the quark mass m, which have been fitted to the pion charge radius and decay constant. 

Inserting the model wave function of Eq. fl6UI) in the LCWF overlap representation of 

GPD in Eq. fISS]) . one obtains 

dK± ^ ^ I \ , / \ + ■ 

X (l-x) Mo{x, k'J Mo{x, k^) ■ 

(65) 



i/"(x,^ = 0,t) = J ij^{x,K'^)ip^{x,K±) 



In the forward limit — > 0, the Melosh rotation matrices combine to the identity matrix 
and one obtains the ordinary parton distribution as momentum density distribution given 
by the square of the momentum-dependent part of the wave function [l7|, i.e. for x > one 

gets 

<x) = J mx,^^)\\ (66) 
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VI. RESULTS AND DISCUSSION 

In this Section the results obtained from the different models described in the previ- 
ous Sections are presented and discussed. Let us first illustrate the actual values of the 
parameters entering the three models. 

For the covariant model (Sec. IIV|) the CQ mass and the pion mass have values m = 220 
MeV and m.,^ = 140 MeV, respectively. It should be pointed out that, for some runs, the 
value = has been used in order to match the vanishing pion mass adopted for the 
Mandelstam-inspired model (see Sec. IV Al) . This change will be adequately emphasized 
whenever applied (in this case the CQ mass is a little bit lowered, i.e. m = 210 MeV). The 
parameter mji present in the pion Bethe-Salpeter amplitudes is fixed through the pion decay 
constant, obtaining = 600 MeV for the sum- form (Eq. fl37j) ) and rriji = 1200 MeV for 
the product-form (Eq. (l38!l ). 

In the Mandelstam-inspired model, as already mentioned, all the parameters are the same 
ones used in j2lj, except for i) wvm = — 1 that yields the relative weight of the instantaneous 



contributions and ii) the CQ mass, m = 200 MeV, i.e. the one adopted in [48| within the 
same approach for the very detailed description of the nucleon em form factors in both the 
spacelike and timelike region. As already mentioned, given the complexity of the calculation 
a simplifying assumption of a vanishing pion mass has been also added. Finally, in the 
VM dressing of the quark-photon vertex (cf Eq. (150|) ) up to 20 isovector mesons have been 
considered in order to have a good convergence even for t = —10 (GeV/c)^. 

In the LFHD model (see Sec. IV Bl) . a CQ mass m = 250 MeV and a wave-function 
parameter (3 = 319.4 MeV have been used in order to reproduce the pion charge radius 
{i^ch = 0.670 ± 0.02 fm) and the pion decay constant 

First of all, the theoretical models have been compared with available experimental data, 
in particular the pion em form factor in the spacelike region. 

In Fig. [3], it is shown the ratio between the spacelike form factors, calculated by using 
our models, and the monopole form factor Fmon = 1/(1 + l^l/'^p) i^p = ■''"''0 GeV). The 
relevance of such a presentation of the form factor is twofold: i) dividing by Fmon allows 
one to avoid the log plot that hinders a detailed analysis, ii) more important, one can 
immediately discriminate between models that produce a divergent charge density at short 
distances and models that do not (cf, e.g., ^]), since their fall-off is more rapid than Fmon- 
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-t (GeV/c) 



FIG. 3: Pion form factor vs — t. Thin dashed hne: covariant symmetric model of Ref. [20(], 
with the momentum dependence of the pion Bethe-Salpeter amphtude given by the sum-form of 
Eq. (j37p . and m-,^ = 140 MeV. Double-dot-dashed line: calculation performed within the LF 



Mandelstam- inspired model (cf Sec. IV Ap . by using an asymptotic pion wave function 



271] with 



0, and adopting a CQ mass of m = 200 MeV (notice that in [21] m = 265 MeV). Thick 



521 ] ■ arbitrarily extended in this figure 



solid line: monopole fit to Lattice data as obtained in Ref. 
from —4 (GeV/c)^ to —10 (GeV/c)^ (see text). Thick dot-dashed line: faster-than-monopole fit 
to Lattice data as obtained in Ref. [52], arbitrarily extended in this figure from —4 (GeV/c)^ to 
—10 (GeV/c)^ (see text). Dot-dashed line: the same as the double-dot-dashed line, but with a 



non perturbative pion wave function, eigenstate of the squared LF mass operator of Ref. |42l ]. 
Dotted line: the same as the thin dashed line, but with the product-form of Eq. (|38p for the pion 
Bethe-Salpeter amplitude. Thick dashed line: LFHD model (cf Sec. IV Bp with a Gaussian pion 
wave function and the proper Melosh rotations. Experimental data: full dots from the collection 



of Ref. 



491 ]: open squares, TJLAB data from Ref. 



For the sake of completeness, we have also displayed two different fits (thick solid and 
dot-dashed lines in Fig. [3]) to the Lattice data as obtained in Ref. j52|. In that paper. Lattice 
data have been extrapolated to the experimental pion mass, and they were described up to 
t = -4 (GeV/c)2 both in terms of i) a monopole function F^"*(t) = 1/[1 - t/M^{ml''y')] 
with MiwPy') = 0.727 GeV and ii) a function with a fall-off faster than the monopole one, 
Le. Fl^^t) = -t/{p M^{mPj'y')]P with p = 1.173 ± 0.058 and M(mP''?'^) = 0.757 ±0.018 
GeV. In Fig. |3] the Lattice results have been arbitrarily extended by using the previous 
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FIG. 4: Pion form factor calculated within the covariant model of Sec. IIVI with and without a 
vanishing pion mass. Solid line: sum- form for the pion Bethe-Salpeter amplitude, Eq. (j37p . and 
rTT-Tr = 0. Dashed line: the same as the solid line, but with = 140 MeV. Dot-dashed line: 
product-form for the pion Bethe-Salpeter amplitude, Eq. ([38]), and = 0. Dotted line: the same 
as the dash-dotted line, but with m-,^ = 140 MeV. Experimental data as in Fig. [3l 

functions from t = —4 (GeV/c)^ to t = —10 (GeV/c)^, with a quite reasonable outcome. 

To show the sensitivity of the covariant model of Sec. UV] upon the change of the pion 
mass, a comparison between calculations performed with a vanishing pion mass and with 
m-TT = 140 MeV is presented in Fig. HI These calculations are helpful in view of the 
following comparisons with the Mandelstam-inspired model, where the value = has 
been adopted. It is interesting to notice from Figs. [3] and H] that the sum-form for the BS 
amplitude is unable to accurately describe the experimental em form factor at high values 
of \t\. 

In order to illustrate the frame dependence of the Fock decomposition of the em form 
factor, in Fig. [5] the valence and non- valence contributions to the pion form factor (Eq. ( 1271) ) 
within the covariant model based on the product-form and = are presented. Such a 
choice for m.^ is suggested (cf Eq. (135|) ) by the need to explore the whole range < |^| < 1. 
The sum of the two contributions becomes ^ independent, and the result is shown in Fig. 
m by the dot-dashed line. Figure [5] allows us to disentangle the valence and non-valence 
contributions. Indeed, different values of ^ correspond to different choices of the frame (let 
us remind that ^ = corresponds to the Drell-Yan frame and = — 1 to the frame where 
A_L = 0). Moreover, it is worth noting that the operator "number of constituents" does not 
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FIG. 5: Left Panel: valence contribution, Fj^" {\^\,t), to the em pion form factor (see Eq. (j25p ). 
evaluated within the covariant symmetric model of Sec. IIVI by using the product-form for the 
momentum-dependent part of the Bethe-Salpeter amplitude (cf Eq. (j38|) ) and choosing = 
for covering the whole range < |^| < 1, according to Eq. (j35p . Right Panel: the same as in the 
Left Panel, but for the non-valence contribution, Fi™\\^\,t), (see Eq. ([26]) ). Note the different 
orientations of the axes in the two Panels, for a straightforward selection of the relevant regions. 

commute with the whole set of the Poincare generators, and therefore a change of frame 
alters the non- valence content. The knowledge of valence and non- valence contributions in 
the plane (^, t) could impose new constraints to models that aim to go beyond the standard 
CQM. 

After completing the analysis of the em form factor within our models, in Fig. [6l the 
isovector GPD for positive x, namely H^^^{x,0,0) = Uv{x)/2 (see Eq. (fT9|) ). is shown as 
a function oi x = Xq (since for ^ = one recovers the longitudinal momentum fraction, 
Eq. ([3])). It should be pointed out that, at this stage of our analysis, no evolution has been 
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applied. The effects of tlie evolution for the parton distribution will be considered elsewhere, 
together with a study of the evolution for the whole GPD. Calculations for the covariant 
model of Sec. [IV] and the LFHD model of Sec. IV Bl are shown in Fig. [6l Notice that the 
Mandelstam-inspired LP model presently allows predictions only for |^| = 1. In order to 
extend to ^ = this approach, a non-vanishing value of and a bare term, besides the 
VMD one, should be considered. Thus one can take into account the contribution depicted 
in Fig. Wia), that produces the valence term in H^^^{x, ^, t), but new, non-trivial parameters 



have to be added fcf the nucleon case in 



48|). 



The comparison in Fig. [6] shows the difficulty of the sum-form, Eq. ( 137|) , for the pion 
Bethe-Salpeter amplitude to give a realistic parton distribution, i.e. to have a vanishing 
value at the end-points. Reminding that, for ^ = 0, the presence of the delta-function in 
Eq. fBOj) and the kinematical relations in Eq. ([1]) impose the correct support [0, 1] for the 
variable x (remind that for ^ = 0, one has x = Xq), the sum-form produces a discontinuity 
at the end-points, i.e. an infinite derivative. It is instructive to correlate such a drawback 
to the one already seen in Fig. [3l where the sum-form is not able to reproduce the em 
form factor at high values of (— t). Indeed, in both cases, the high momentum part of the 
valence component of the pion state is involved. As a matter of fact, for a; = and x = 1, 
the intrinsic three-momentum becomes infinite (cf Eqs. (16T!) and (!62l) ). and therefore small 
distances are involved, just as in the case of the tail of the em form factor, where the infiuence 
upon the small-r part of the pion wave function is felt. The more realistic behavior of the 
product-form ( 138|) . can be ascribed to a \k±\ fall-off like the one dictated by a BS kernel 
dominated by a one-gluon-exchange, as already pointed out in Sec. IIVI An important, final 
remark is the clear shift towards small x of the curves evaluated within the covariant model, 
while the prediction obtained within the LFHD model is symmetric with respect to x = 1/2. 
Such an interesting difference could be explained by the fact that the full covariance of the 
model of Sec. HV] together with its dynamical content, related to the adjusted parameter 
rriR, could take into account some effects beyond the pure qq component of the pion state. 
First, one should note that the valence component, Eq. ( l39i) . generates a quark distribution 

y less than 1: for the sum-form 
20||). Then, by using the Fock 
decomposition of the pion state (see, e.g., {4] for a general discussion), one immediately 
recognizes contributions from both the qq component, (i.e. the valence component) and 



symmetric with respect to x = 1/2 and a probability definite^ 
Pvai = 0.78 and for the product-form Pyai = 0.84 (see also 
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FIG. 6: Isovector GPD H^^^{x,0,0), equal to half of the parton distribution (see Eq. (|19p ). vs x. 
Thin dashed line: covariant model of Sec. IIVI calculated by using the sum- form, Eq. (j37p . for the 
pion Bethe-Salpeter amplitude and m.,^ = 140 MeV. Dotted line: the same as the thin dashed line, 
but for the product form, Eq. (j38|) . Thick dashed line: LFHD model of Sec. IV Bl with a Gaussian 
pion wave function and the proper Melosh rotations. The variable x, given in Eq. ([T|), coincides 
with the usual LF longitudinal fraction Xq, since ^ = (see text below Eq. (jlOp ). 



from other components with more constituents (see e.g. the instantaneous contributions in 
Eqs. (IBIOI) . (1B11I) and ( 1B12I) and the analysis in 55|]). Thus, the active quark shares the 
longitudinal momentum of the pion with more than one spectator parton, belonging to the 
Fock space configuration beyond the valence one. Therefore, the shift toward values of x 
less than 1/2 is expected, since our covariant model contains more physical effects than the 
basic one. In particular, for a non vanishing pion mass the average longitudinal momentum 
fraction for the sum-form is < Xq >~ 0.483 and for the product-form is < Xq >~ 0.471, i.e. 
quite similar, but a little bit different from 1/2. As a simple cross-check we have reobtained 
those values also from v42^°(0) =< Xg > (cf Eq. ([22]) with j = and Eq. ffTSl) ). 

A more detailed analysis of the parton distribution can be achieved by using the 
chiral-even TMD distribution, fi{x, |k^|), see Eq. fl20l) . In Fig. [71 the TMD distribu- 
tions calculated within the covariant model by using the different BS amplitudes of Eqs. 
( 1371) and ( l38l) are shown. In order to avoid log plot, /i(x, |kj_|) has been divided by 



G{\k±\) = 1/(1 + |fc^|7m2)4, (with rup = 770 MeV)). Clearly, the product-form has a |k_L| 
fall-off faster than the sum-form does, i.e. low transverse-momentum partons are favored in 
the first case. 
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FIG. 7: Transverse-momentum dependent function, |fc_Lp)/G'(|fc_L|), with G(|fc_L|) = 1/(1 + 
\k / m?p)'^ . Left Panel: sum-form of the Bethe-Salpeter amphtude (see Eq. (j37p ). Right Panel: 
the same as in the Left Panel, but for the product-form of the Bethe-Salpeter amplitude (see Eq. 
(j38|) ). The normalization is given by dx J d'k± fi{x, |k_Lp) = 1, and kperp means |k^|. 

The analysis of both the generalized form factors involved in the second moment of the 
isoscalar pion GPD, i.e. o°(^) and ^2l°(^) (cf Eq. (!22|) with j = 0), has to be performed 
necessarily within the covariant analytic model of Sec. IIVI This is obvious if we look at Eq. 
([17]), where the polinomiality imposes a square dependence upon ^, and therefore one needs 
a model that covers an extended range for the variable ^. Indeed, for each value of t, we 
have first numerically checked the parabolic behavior against ^, and then we have extracted 
the coefficients of the parabolic fit getting the values of o°(^) and ^2l°(^)- Figure E] shows 
a comparison between i) recent results from Lattice QCD, extrapolated to the physical pion 
mass 0, [5^ , ii) our covariant calculations evaluated with both = and = niphys by 
using the sum- and the product-form for the BS amphtude (Eqs. (|37|) and fl38|) ) and iii) the 
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LFHD result (see Sec. IV Bp for Al'^^it) only, since this approach at the present stage allows 
one to perform calculations exclusively for = 0. Indeed, the ratios o°(^)Mo o°(0) and 
^2^^i't)/^2^^{^) are presented in order to get rid of the evolution (see Ref. [13| for a detailed 
discussion of this issue). The Lattice calculations are described through a monopole form, 
1/(1 — t/M|j), as obtained in 53| from the analysis of their Lattice data, without evolution 



and with evolution in the MS scheme at the scale fx = 2 GeV. In particular, we have used the 
following values: Mg.o = 1.329 ±0.058 GeV and M2,2 = 0.89 ±0.25 GeV, corresponding to an 
analysis of the Lattice data that satisfies the low-energy theorem, i.e. A^=0(0) = -4A^=°(0). 
The uncertainties on the previous masses generate the shaded areas in the Left and Right 
panels in Fig. [HI 

Unfortunately, i) the available range of {—t) (we refrained to enlarge the interval as we did 
in the case of the em form factor, since we do not have experimental data yielding confidence 
in an arbitrary extension of the monopole fit) and ii) the large uncertainties in the Lattice 
calculations of ^2^2° do not allow us to elaborate too much on the comparison between our 
phenomenological models and the Lattice results. On the other hand, for large values of |t| 
the calculations obtained by using the covariant model with the product-form and m^r = 140 
MeV could give some insight on the expected behavior of the Lattice calculations, since one 
could argue that the covariant model with the product-form phenomenologically contains 
at some extent dynamical features typical of QCD, like the one-gluon-exchange dominance 
at small distances. In order to complete the information, in Table I the values of A2"^^{0) 
and ^2'2°(0) shown. It is worth noting that while the Lattice calculations largely fulfill 
the low-energy theorem, as already mentioned, our calculations do not. Furthermore, it 
should be pointed out that for small t the disagreement between Lattice data and the calcu- 
lation with the covariant approach at some extent is an expected one, since the mechanism 
responsible for the confinement is not present in our model, and therefore we have a free 
propagation of the qq pair. A possible solution could be elaborated following the suggestion 
in Ref. 56[ , where a covariant model without the disturbing free propagation of the qq pair 
was proposed and applied to the em decays of the vector mesons. 

The previous figures have illustrated "integral" properties of the pion GPD's, like em form 
factor and the generalized ones, or the parton distribution, i.e. H^^^{x, 0, 0). In the following 
figures, the isoscalar and isovector GPD's are shown in the plane (x, t) with — 1 < x < 1 and 
— 10 (GeV/c)^ ^t<0, but with fixed values for ^, as dictated by the two phenomenological 
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TABLE I: Gravitational form factors at t = 0, (cf Eq. (j22p with j = 0), obtained i) within the 
covariant model of Sec. IIVI and both the sum- and the product-form for the BS amplitude, Eqs. 
(I37p and (I38p : and ii) from the Lattice data of Ref. 





Sum 

m,r = 


Sum 

TTItt — ^phys 


Product 
m,r = 


Product 

TUty — TTT'phys 


Latt. no evol. 


Latt. with MS evol. 

TUty — TTT'phys 




0.4828 


0.4833 


0.4707 


0.4710 


0.365 


0.261 




-0.0307 


-0.0272 


-0.0357 


-0.0327 


-0.092 


-0.066 




o 0.8 -\ 



< 0.2 - 




-t (GeV/c) 



FIG. 8: Left Panel: the ratio ^2^0° (0/^2^0° (0) ' involving the generahzed form factor ^2^0° (*) ^^^^ 
appears in the second moment of the isovector GPD H^^^ (cf Eqs. (|17p and (j22p l as a function of 
t. Solid line: sum- form for the pion Bethe-Salpeter amplitude, Eq. (j37p . and m,r = 0. Dashed line: 
the same as the solid line, but with m,r = 140 MeV. Dot-dashed line: product-form for the pion 
Bethe-Salpeter amplitude, Eq. (j38p . and rn-jr = 0. Dotted line: the same as the dash-dotted line, 
but with TTT-Tr = 140 MeV. Thick long-dashed line: LFHD model (cf Sec. IV Bp with a Gaussian pion 



53] (see 



wave function and the proper Melosh rotations. Shaded area: results from Lattice QCD 
text). Right Panel: the same as the Left Panel, but for A^^f (t) / A^^f (0) . 



models, namely |.^| = 1 for the Mandelstam-inspired model (Sec. |V]) and ^ = for the 
LFHD model (Sec. IV Bp . respectively. The covariant model (Sec. lIVp . in its two versions 
for the momentum dependence (Eqs. fl37|) and (l38l) ). will be compared to the results for 
the two phenomenological models, that, in some sense, represent two extrema, in the Fock 
language: the first model is basically related to the non-valence (ERBL) region, the second 
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one is related to the valence (DGLAP) domain. In order to cover the whole range of ^ for the 
given interval of t (i.e. —10 (GeV/c)^ ^ t < 0) the covariant model has been evaluated by 
assuming = 0, as already pointed out. Finally, let us stress that the GPD's are divided 
by Fmon, as in the case of the em form factor, for avoiding log plot and for emphasizing 
as many details as possible. In Fig. [9], the results of the covariant symmetric model are 
shown for |^| = 1, in order to be compared with the calculations performed by using the 
Mandelstam-inspired model, presented in Fig. [TOl We remind that the phenomenological 
model has a photon-quark vertex dressed by a microscopical VMD, as discussed in Sec. |Vl 
In Fig. [TT], the no-helicity flip GPD's of the covariant symmetric model are shown for ^ = 0, 
allowing a comparison with the calculations performed by using the LFHD model, presented 
in Fig. [121 For ^ = 0, where only the valence component is acting and | x| — "^q^ a nice 
feature, stemming from the figures of both isoscalar and isovector GPD's, is shared by all 
the presented models: in the limit of large \t\ the collinearity clearly emerges, as shown by 
the migration of the maximum (minimum) value from |x| ^ 1/2 for t = toward |x| ~ 1 for 
\t\ oo. Such a behavior can be easily understood in the LFHD model, since non vanishing 
contributions to GPD (cf Eq. ( l65l) ) can be obtained if k'j_ in Eq. ( I59|) does not depend too 
much from A_|_, namely a; ~ 1 (notice that for x exactly 1, the free mass blows up and the 
wave functions become vanishing, as well as GPD's). Correspondingly, for |^| = 1, where 
only the non- valence component is acting, the relevance of the x region around ±1 can be 
explained by the pair-production mechanism. For simplicity, let us consider large values 
of that amount to large values of = (remind that in the Breit frame A*^ = 0). 
Then, using A"*" = + = kzq + kzq ~ 2^^^ > (given our choice for the sign of A+), 
and the fact that each quark in the pair is almost on its mass-shell, we can approximate 



2k~^ = kg — kq 2Eg = 2\J m? + Thus, one can see that, when ^ = —1, x becomes 

close to 1 for A^ >> m, since x = k^/P^ = 2k^/A'^ ~ Eq/k^q — * 1. The case x = —1 can 
be obtained for A+ < 0. 

Finally in Fig. [131 the isovector GPD, evaluated within the covariant model adopting 
the product-form of Eq. ( 138|) and m^r = 0, is shown for the case x = \C,\ and > t > —10 
(GeV/c)^. This kinematical region, where the transition from DGLAP to ERBL regimes 
occurs, should be relevant for the experimental studies of the single spin asymmetry (see, 
e.g. the discussion in ^\). Let us notice that in our covariant analytic model the GPD is 
continuous at a; = |,^|. 
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H^i=o(x,l^l=l,t)/F^Jt) H^i=i(x,l^l=l,t)/F^Jt) 




FIG. 9: Upper Left Panel: Isoscalar no-helicity flip GPD from the covariant symmetric model of 
Sec. Ml with the sum-form for the Bethe-Salpeter amplitude (Eq. ()37p ) at |^| = 1. The value 
of ^ is fixed by using = (cf Eq. (|35p ). for the sake of comparison with the microscopic 
model of Sec. |Vl whose results are shown in Fig. [TOl On the z-axis the ratio with respect to 
Fmon = 1/(1 + l^l/''^p) is presented. Upper Right Panel: the same as in the Upper Left Panel, but 
for the isovector GPD. Lower Panels: the same as in Upper Panels, but for the product-form for 
the Bethe-Salpeter amplitude (see Eq. (1381)'). 
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H^i=o(x,l^l=l,t)/F^Jt) _ H,i=i(x,l^l=l,t)/F^Jt) 




FIG. 10: Left Panel: Isoscalar no-helicity flip GPD from the Mandelstam-inspired model of Sec. 
|V]at 1^1 = 1 (see text). Right Panel: the same as in the Left Panel, but for the isovector GPD. 

From the 3D plots, one can see that the covariant model, in the version with the product- 
form for the momentum-dependent part of the BS amplitude, is able to reproduce quite sat- 
isfactorily the GPD's evaluated within the two phenomenological models, the Mandelstam- 
inspired and the LFHD ones, and therefore one could argue that it contains the main in- 
gredients for a realistic descriptions of the constituents inside the pion. In view of this, it 



appears challenging to test the covariant model (or its refinements 57| based on the Nakan 



ishi representation, see, e.g. |36|, for a recent applications to a bosonic system) of the BS 
amplitude, in comparisons with experimental data, whose analysis requires the knowledge 
of the pion GPD's. 
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FIG. 11: Upper Left Panel: Isoscalar no-helicity flip GPD from the covariant symmetric model of 
Sec. llVt with the sum-form for the Bethe-Salpeter amplitude (Eq. (f37|) ) at ^ = 0, and m,r = 0. 
Upper Right Panel: the same as in the Left Panel, but for the isovector GPD. Lower Panels: the 
same as in the Upper Panels, but for the product-form for the Bethe-Salpeter amplitude (Eq. (|38]l ). 
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H i=o(x,£=0,t)/F (t) H i=i(x,£=0,t)/F (t) 




FIG. 12: Left Panel: Isoscalar no-helicity flip GPD from the LFHD model of Sec. |VB] at ^ = 
(see text). Right Panel: the same as in the Left Panel, but for the isovector GPD. 

VII. CONCLUSION 

In this paper, we have investigated the no-helicity flip Generalized Parton Distributions 
of the pion by using three models, based on a description of the pion where constituent 
quarks with masses between 200 MeV and 250 MeV are considered. In particular, we have 
evaluated the isoscalar and isovector GPD's adopting a covariant, analytic model and two 
Light-front phenomenological models. It is important to notice that the first model, based on 
4D Ansatzes for the Bethe-Salpeter amplitudes, allows us to explore the whole kinematical 
domain of the three variables x, ^ and t upon which the GPD's depend, while the others 
two are, presently, constrained to a given value of ^. The second model, the Mandelstam- 
inspired model of Sec. |Vl is a natural extension of the approach proposed in Ref. {21] for 
a successful investigation of the em form factor of the pion in both the space- and timelike 
regions. Main features of the model are: i) a microscopical Vector Meson Model dressing 
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H^i=i(x,l^l=x,t)/F^Jt) 




X 



FIG. 13: Isovector no-helicity flip GPD from the covariant symmetric model of Sec. IIV( with the 
product-form for the Bethe-Salpeter amphtude (Eq. (f38|) ) at |,^| = x, and m-jr = 0. 

for the quark-photon vertex and ii) proper Ansatzes for the 3D LP projection of the BS 
amphtudes of both pion and vector mesons, taken as the eigenfunctions of a LP square mass 



operator [42|]. As in [2l|, the assumption = is added, and this simphfication allows 
calculations of the GPD's only for the value |.^| = 1 (cf Eq. fl5^ ). namely the non-valence 
region covers the whole range 1 > a; > — 1. 

On the contrary, the LPHD model of Sec. IVBl is based on a Poincare covariant description 
of the pion, with a proper treatment of the spin wave functions, due to the presence of the 
Melosh rotations. The momentum part of the pion wave function is given by a Gaussian 
function, that contains the dynamical input of the model through two adjusted parameters. 
A bare quark-photon vertex is assumed. It is worth noting that the model yields a description 
of the GPD's for ^ = 0, i.e. the valence region can be investigated. 

The covariant symmetric model of Sec. llVt based on a Mandelstam formula for matrix 
elements of the operators yielding the isoscalar and isovector GPD's, allows us to have close 
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expressions for the physical quantities, since analytic forms for the momentum-dependent 
part of the Bethe-Salpeter amplitude are adopted and a bare quark-photon vertex is assumed 
as well. Such a covariant model can be applied for any value of x, ^ and t, and therefore can 
be used for interpolating between the two phenomenological models. A peculiar feature is 
given by the presence of instantaneous terms, both in the valence and non-valence regions, 
since we fully take into account the analytic structure of the BS amplitude. 

The comparison with the em form factor (Fig. [3]) suggests that the covariant model with a 
sum-form of the BS amplitude has a non realistic increasing behaviour with respect to Fmon, 
for large which leads to a divergent density at short distances, while the version with 
a product-form together with the LFHD model decrease more rapidly than Fmon- Finally, 
the Mandelstam- inspired model and the Lattice results (red curve in Fig. [3]), arbitrarily 
extended from t = —4 (GeV/c)^ to t = —10 (GeV/c)^, given the analytic form proposed in 



52l | for extrapolating the Lattice data to the physical m^, show a moderate decreasing with 



respect to Fmon, for large \t\. Such a comparison for the em form factor and the analysis of 
the parton distribution in Fig. [6] point to the relevance of the behavior of the pion valence 
function (or better the momentum part of the BS amplitude) for large transverse momentum. 
In particular the product-form, that has a behavior at large transverse momentum |k i | 
compatible with the one suggested by the one-gluon-exchange dominance (see, e.g. |36|). 
seems to give a consistent description of both the tail of the em form factor and the end-point 
fall-off of the parton distribution. With respect to this finding, more details can be gained 
from the investigation of the chiral-even transverse-momentum dependent distribution, as 
shown in Fig. [3 

Another important step in the characterization of the covariant model is given by the 
comparison of the generalized form factors with the Lattice results. For the present, the 
comparison is restricted to the gravitational form factors, A2"^^{t) and A2'^^{t), that appear 
in the second moment of the isovector GPD, H^^^, (cf Eqs. ( fTTl) and f l22l) ). Indeed, for 
^2 o°(t) we have presented results from both our covariant model and the LFHD approach, 
while for A2'^^{t) only the covariant calculations are available (let us remind that calculations 
with ^ 7^ are necessary for disentangling both form factors). Unfortunately, since Lattice 
data have been obtained in a t-interval not too wide and are affected by large uncertainties, 
one cannot yet draw stringent conclusions from the comparison shown in Fig. [HI However, 
the encouraging agreement between model calculations and Lattice data for both ratios. 



38 



A^=°(t)/A5=°(0) and A^=°(t)/A^=°(0), suggests to extend our anal ysis also to the spin-flip 
GPD's, since Lattice results are available for the lowest moments 5J], in order to explore 
the onset of the dominance of a one-gluon-exchange mechanism for a light hadron. 

To complete our analysis, we have studied the GPD's in the {x, t) plane for flxed values 
of ^, i.e. 1^1 = 0,1, X. These values are representative of different, interesting cases. The 
flrst one, ^ = 0, involves contributions to GPD's only in the valence region, while the 
second one involves contributions only from the non-valence one. Finally, the case |^| = x 
illustrates the transition from the DGLAP region to the ERBL one. The covariant model 
can explore the whole 3D space of the variables (x, ^, t) and it is compared with LFHD 
model for ^ = 0, and with the Mandelstam-inspired model for |^| = 1, while for |,^| = x 
shows a smooth transition from DGLAP region to the ERBL one, given the continuity of 
the model. It should be pointed out that the covariant model with the product-form for 
the Bethe-Salpeter amplitude exhibits an overall agreement with the Mandelstam-inspired 
model, for |^| = 1, and with the LFHD model, for ^ = 0. Therefore, from these flndings 
one could conjecture that the general shape, illustrated by the previous covariant model 
and the phenomenological ones, is a typical feature of the pion GPD's, dictated from both 
kinematical arguments (cf the discussion at the end of Sec. IV1|) and the dynamical input 
reflected by the proper fall-off of the momentum distribution (cf the one-gluon exchange 
dominance at short distances). 

Further analyses, to make more and more realistic the models presented in this paper, 
are in progress. 
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APPENDIX A: KINEMATICS 

Following the notations of Fig. [H where A"*" > 0, one obtains from Eq. ([1]) that > ^. 

In the valence region, for a quark, one has: i) in the initial state, > — > 0, 

i.e. > > (notice that necessarily the spectator constituent is an antiquark, 

since > A;^ — P^) and then 1 > x > —C,', ii) in the final state, p'^ > k^ + > 0, i.e. 

P+ > k^ > — A+/2, and then 1 > a; > ^. Therefore in the valence region, one gets the 
interval 1 > x > — ^, and given our choice for ^ one has 1 > x > |^|. 

For an antiquark in the initial pion, the four- momentum is k + A/2, while the spectator 
quark has four-momentum k + P. In the final pion, the antiquark four-momentum is A; — 
A/2. The antiquark plus components are negative both in the initial and in the final 
pion. Therefore i) p^ > —{k^ + A+/2) > that leads to > x > — 1 and ii) p'^ > 
— {k^ — A+/2) > 0, i.e. — ^ > x > —1. Summarizing, for an antiquark in the valence region 
one finds — > x > —1. 

In the non-valence region, one has to deal with a qq production, i.e. > fc"*" — /2 
and A;+ + A+/2 > (see Fig. [21), and those constraints translate into ^ < x < — ^. The qq 
annihilation is prevented by the choice of a positive A+. In order to have general extrema, 
holding for both positive and negative A+, one can write |^| > x > — 

APPENDIX B: INTEGRATION ON k' 

In this Appendix, the no-helicity flip GPD for the symmetric covariant models (see Sec. 
irv]) calculated using Eq. (HO!) and the momentum dependent part of the BS amplitude, 
given by Eqs. (|37D or (|38l). 

The evaluation of the trace in Eq. (jUJ) can be simplified according to the decomposition 
of the Dirac propagator shown in Eq. ( H2l) and reminding that [7"*"]^ = 0. By introducing 
the variable k = P — k, one has 

Tr[0^{K')] = Tr {(^ + m) (/ - ^ + m) 7+ (;^ - ^ + m)} = 

= Tr[O^K-J] + Tr {7+ [(/ - + m] 7^ [{^ - + m]} = 

-{P+-K+) [ip'~KU-Kon-m^]}+A {k--k-^;) {p'+-k+) (p+ - /«+) (Bl) 
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where 

Tr[0+{^-n)] = Tr + m) [(/ - + m] 7+ - ^)on + ^]} (B2) 

After performing the scalar products, one gets 

Tr[0+{K-)] = 4 p'+p+ - |AxP + 2A+ ■ + 

+4 {k- - - (p+ - K^) (B3) 

Given the simple expression adopted for the momentum dependence of the BS amplitude 
(see Eqs. f l371) and fl38l) ). the analytic integration on fc~ can be easily performed in Eq. PUI) . 
By using the LF variables (i.e. (i^K dK^dn~dn^/2) one obtains 



i^n^,e,^) = -^iV.7^ / ^^^^^^^ ^ [P^l -a:) - 



1 1 1 



X 



AiK,p') AiK,p) (B4) 
where 

^ ~ '^)on = 7ZZ "In [V — 1^ ~ 



K+ (p+-fi:+) (p'+-/€+) 

(B5) 

In the integration over the minus component k,^ one faces with the following six poles 
(coming from the BS amplitudes and the Dirac propagators) 



^1(2) -^on{R) \^ 



\p^-K+) 



S(6) =P' ~ (P - '^)on(R) ' " ^^1+ _ 



where k^, (p — and (p' — k)^ can be obtained from the corresponding quantities in Eq. 
(IB5I) by substituting m — > niR. Notice that K2 can appear both as a single and as a double 
pole. 

It is easily seen that the analytic integral flB4p is not vanishing only if > k+ > 0. 
Furthermore we can recognize two subinterval i) p+ > k+ > 0, or valence region, and ii) 
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p'"*" > > p"*", the non- valence region. Let us stress that Eq. (lB4p is vanishing for x < 
since in this case k+ = P+ {1 — x) > P+ (1 + |^|) = p'^- 

In the valence region, only the poles and k,2 belong to the lower semiplane. In the 
non-valence region, only /tg and Kg belong to the upper semiplane. 

To obtain the no-helicity flip GPD in the valence region, let us integrate over closing 
the contour in the lower semiplane. The contribution from reads as follows 



6[P+{l-x)-K+] 



A(«:,P')L- 



[P - l^on -{P- f^)on] [P' - f^on " (P' " f^)on] 

where for the sum-form, Eq. (^7} . one has 

H^^p)L. = Ci • ^ 



(B7) 



+ 



1 



{p+ - K+) [p- - K-„ -{p- k)^] 

and for the product-form, Eq. fl38l) . one has 

1 1 



fB8) 



{p+ - K+) [p - K„„ -{p- fi^if^on - ^^r) 



(B9) 



For the sum-form, the pole Kj^ generates a contribution as a single pole and a contribution 

ClS db double pole. 

The single-pole contribution is given by 



4(27r)^ 



dn 



k'^R Ko„) [p 



6[P+{l-x) - K+] 

{K+y{p+ - K+){p'+ - K+) 
X 



Tr[0+(«:^)] X 



(P - '^)on)] b' -f^R- iP' - ^^)on\ 



1 



{p+ - K+) - - (p - k)^] - /€+) [p'- -Kj^- ip' - K,)j^] 

and the double-pole contribution is given by 



(BIO) 



W^v-'^'-'/ 4(2vr)3 
d {Tr[0+{K-)] 



Cf / c/kj 



S[P+{1-X) - K+] 
{K+f{p+ - K+){p'+ - K+) 



dtv 



{l^ - l^on) [P -{P- K.)on)] [P' - - (P' - l^)on] 



(Bll) 



For the product-form, the pole Kn generates only a double-pole contribution, given by 



4(27r) 



3 ^2 



dn^- 



6[P+{l-x) - /^^ 



K+)3{p+ - K+)^{p'+ - K+)2 



X 
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d f Tr[0+{n-)] 1 ^ 

I (l^^ - l^on) [P' - - (P - /^)on)] b'" -K- - {p' - k)-J 



[p'-K -(p-k)^)] [p'--K -{p'-K)f^] 



(B12) 



The contribution in the non-valence region can be evaluated by considering the poles 
and Kg . In particular the contribution from Kg has the same form for both choices of the 
BS amplitudes, i.e. 



[p'- - (p' - k)-„ - K-„] [p- - p'- + (p' - k)-„ - - k)^„] 
while the contributions from Kg , reads differently for the sum-form, viz 

rru ^ . ^ 7^ ^ /• , , , 5 [P+(l - X) - K+] 

^^^'^ 4(27r)3 J _ _ 

A(/t,p)|p,-_(p,_^)- 

[p- - p'~ + {j>' - n)-^-{p- /t)-„] 
and for the product-form, viz 



X 



(B13) 



V2 ^ 



(B14) 



rru . . N ^ f , r ,+ S[P+{1-X)-K^ 



Xt -— T T ^ — — - -1 (B15) 

- + (p' - k)j^ ~{p- k)„„J [(p' - {p' - k)^ - 

Summarizing, for the sum-form one has 

H^ix, t) = eix - 1^1) ^(1 - x) (x, e, t) + t) + Hl^^^.ix, e, t)] + 

-x) em+x) +i/r„.)6(x,e,t)] (B16) 

with HJ^nv)6 given by Eq. (]B14p . while for the product-form one gets 

H-{x, e, t) = 9{x - lei) ^^(1 - x) [Hi^,^{x, e, t) + t)] + 

- a:) Om+x) [i/("„,)5(x,e,t) + i^r„.)6'(x,e,t)] (B17) 
with i/«^)g, given by Eq. 
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APPENDIX C: ELECTROMAGNETIC FORM FACTOR 



The pion electromagnetic form factor is defined by 



1 



2 P^ 



(7r+(p')|J(0)-n|7r+(p)) 



dx H^-^{x,i,t) 



dx 



dz 



(CI) 



2=0 



where the range of x has been extended from [—1,1] to [—00,00], since H^^^{x,C,,t) is 
vanishing outside the support [—1, 1], given the presence of the delta-function in Eq. (140|) 
and the kinematical relations in Eq. ([1]) (see, e.g. [l, 4|). 
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